m 


PASSIVE  INFRARED  RANGING 


THESIS 

Norman  K.  Leonpacher 
Captain,  USAF 

AFIT/GEP/PH/83D-5 


!  This  document  has  been  cpp’cvcd  j  • 

j  for  public  release  and  :sal*.:  it?  I 
i  distribution  is  unlimited.  ! 

i  v  r 

DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


7*  ^LECTB. 

,  .'EQ  2  2  1S84 


Wright-Patterson  Air  Force  Base,  Ohio 


84  02  21  175 


AFIT/GEP/PH/83D-5 


PASSIVE  INFRARED  RANGING 

THESIS 

Norman  K.  Leonpacher 
Captain,  USAF 

AFIT/GEP/PH/83D-5 


!  Thi3  document  has  been  approved 
i  for  public  release  and  sole.  it 
distribution  is  imbn-br  i 


OTIC 

f^ZLECTE| 

FEB  2  2  1984 


A 


AFIT/GEP/PH/83D-5 


Preface 


The  purpose  of  this  study  was  to  show  that  the  range  of 
an  infrared  source  can  be  estimated  by  analyzing  its  spec¬ 
tral  signature.  Although  a  limited  number  of  source  spectra 
were  available  to  develop  and  test  the  ranging  algorithm,  a 
method  has  been  established  that  can  be  expanded  upon  if 
more  data  becomes  available. 

The  VAX  computer  system  at  Air  Force  Wright  Aero¬ 
nautical  Laboratories  (AFWAL) ,  Wright-Patterson  AFB,  Ohio, 
was  used  extensively  in  developing  the  algorithm.  Spectral 
data  files,  and  computer  programs  to  read  these  files,  were 
available  before  this  study  began.  Many  new  programs  that 
have  been  written  are  extensions  of  these  existing  programs. 
Others  are  straightforward  applications  of  the  ideas  in  the 
text,  and  none  are  significant  enough  to  be  included  in  the 
report.  File  copies,  however,  are  available  through  Dr.  R. 
B.  Sanderson,  AFWAL/AAWP-1 ,  Wright-Patterson  AFB,  Ohio. 

Studies  of  this  magnitude  would  not  be  possible  in  the 
limited  time  available  without  the  assistance  of  others.  I 
would  like  to  thank  Dr.  R.  B.  Sanderson  of  the  Air  Force 
Wright  Aeronautics  Lab  for  his  guidance  and  insight  that 
were  instrumental  to  the  success  of  this  project.  I  would 
also  like  to  thank  Nick  Pequinot  for  his  assistance  with  the 
AFWAL  computer.  A  word  of  thanks  is  also  owed  to  my  faculty 
advisor,  Maj  Dave  Graham,  for  his  help  throughout  the 
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project.  And  finally,  I  would  like  to  express  deep  grati¬ 
tude  to  my  wife  Carolyn  for  her  patience  and  inspiration 
that  are  indeed  an  integral  part  of  this  thesis. 


Norman  K.  Leonpacher 
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\  7'  Abstract 

V  I 

The  tfange  of  an  infrared  source  was  estimated  by  ana¬ 
lyzing  the  atmospheric  absorption  by  COQ  in  several  wave- 
•  ^ 

length  intervals  of  its  spectrum.  These  bandpasses  were 
located,1  at  the  edge  of  the  C03  absorr'  ion  band  near 


.  J  fV  |  v 

2300  c^*  (4.3^)* 


A  specific  algorithm  to  predict  rang*  was  determined 


based  on  numerous  computer  generated  spectr 


When  tested 


with  these  spectra,  range  estimates  within  0.8  km  were 
obtained  for  ranges  between  0  and  18  km. 

Accuracy  decreased  when  actual  source  spectra  were 
tested.  Although  actual  spectra  were  available  only  for 
ranges  to  5  km,  63%  of  these  spectra  resulted  in  range 
estimates  that  were  within  1.6  km  of  the  actual  range. 
Specific  spectral  conditions  that  affected  the  range  predic¬ 
tions  were  found.  Methods  to  correct  the  deficiencies  were 
discussed. 

Errors  from  atmospheric  variations,  and  the  effects  of 
background  noise,  were  also  investigated.  Limits  on  accu¬ 
racy  and  range  resolution  were  determined. 


I.  Introduction 


Background 


Any  aircraft  that  flies  into  a  combat  environment  is 


likely  to  become  a  target  for  enemy  weapons.  Because  of  its 


lethal  warhead  and  precision  guidance,  the  surface-to-air 


missile  is  probably  the  most  deadly  weapon.  When  the  rocket 


motor  of  the  missile  burns,  it  emits  intense  radiation  in 


the  infrared  region  of  the  electromagnetic  spectrum.  Infra¬ 


red  receivers  are  being  developed  to  provide  reliable  warn¬ 


ing  of  a  missile  launch. 


In  addition  to  providing  a  launch  warning,  it  may  be 
possible  for  detectors  to  give  an  indication  of  how  far  away 


an  infrared  source  is.  Radiation  is  attenuated  as  it  trav¬ 


els  through  the  atmosphere,  and  its  spectral  profile  is 


altered.  The  constituents  of  the  atmosphere  absorb  at  dif¬ 


ferent  frequencies,  and  the  amount  of  absorption  depends  on 


transmission  range  and  the  concentration  of  each  particle. 


If  a  particular  atmospheric  model  is  assumed,  the  spectral 


signature  of  a  detected  infrared  signal  can  be  compared  with 


a  known  spectral  signature  to  give  an  indication  of  trans¬ 


mission  range. 


Problem 


Specific  frequency  bands  to  use  for  comparison  of 
infrared  signatures  are  sought.  When  two  or  more  bands  are 
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compared,  statistical  distributions  are  formed  that  repre¬ 
sent  spectra  of  sources  of  the  same  range.  The  separation 
between  these  distributions  is  optimized  by  adjusting  the 
frequency  and  bandwidth  of  each  spectral  band.  A  ranging 
algorithm  that  predicts  the  range  of  an  unknown  source  is 
determined  based  on  existing  target  spectra  and  an  atmos¬ 
pheric  transmission  model  called  LOWTRAN  (Ref  1). 

Approach 

In  order  to  develop  the  computer  programs  required  to 
analyze  spectral  signatures,  targets  are  initially  assumed 
to  be  blackbody  radiators.  The  radiant  exitance  for 
blackbodies  is  easily  calculated  from  Planck's  law  (Ref 
2:126),  and  is  spectrally  multiplied  by  LOWTRAN  trans¬ 
mittance  to  simulate  range.  The  resulting  spectra  allow 
development  of  the  analysis  under  controlled  input  con¬ 
ditions  . 

An  example  of  the  radiance  of  a  1000  K  blackbody 
between  2000  and  2500  cm-^  (4-5  ftm )  multiplied  by  two 
different  LOWTRAN  ranges  is  shown  in  Figure  1.  The  atten¬ 
uation  of  the  signal  between  2200  and  2300  cm-^  is  from  CO2 
in  the  atmosphere.  The  total  amount  of  attenuation  is 
related  to  the  range  of  the  source,  but  so  is  the  relative 
attenuation  between  frequencies.  An  indication  of  range  can 
be  obtained  by  comparing  the  relative  signal  intensity 
between  2100  and  2200  cm”l  to  the  intensity  between  2200  and 
2300  cm-1.  If  the  intensities  are  nearly  equal,  the 


Figure  1.  Intensity  of  1000  K  Blackbody  at  0,  1,  and 
10  km 


radiation  could  not  have  traversed  far  through  the  atmos¬ 
phere  so  the  target  is  at  a  close  range.  If  the  intensity 
of  the  higher  band  is  small  compared  to  the  intensity  of  the 
lower  band,  the  source  must  be  at  a  long  range. 

To  compare  several  source  spectra,  a  specific  spectrum 
is  represented  as  a  point  on  a  set  of  coordinate  axes  whose 
values  are  the  intensities  in  each  frequency  band.  Other 
spectra  are  likewise  plotted,  and  if  all  the  spectra  are 
similar  then  the  points  should  cluster.  Distributions  are 
formed  by  signals  that  represent  sources  of  the  same  range. 

The  frequencies  used  in  the  band  comparisons  are 

-  3  - 


/*•*.’ -\v„V * *  . «v#*.  * . ** .  •*.  -  .  •*.  ■*. **.  •**  -  * «/ %•  ■  •  «_*i \*V * 


adjusted  to  provide  maximum  separation  between  each  range 
distribution.  Once  optimum  separation  is  obtained,  ranging 
algorithms  are  determined  based  on  the  statistical  quan¬ 
tities  that  define  the  distributions. 

Assumptions  and  Limitations 

LOWTRAN  is  a  computer  code  that  predicts  the  trans¬ 
mittance  of  the  atmosphere.  To  obtain  a  prediction,  a 
specific  atmospheric  model  must  be  assumed.  The  path 
profile,  including  source  and  sensor  altitudes,  must  also  be 
specified.  In  this  analysis,  only  three  specific  models 
representing  mid  latitude  summer,  mid  latitude  winter,  and 
the  1962  U.S.  Standard  were  used.  In  addition,  the  path 
profile  was  restricted  to  sea  level  altitudes  for  both 
source  and  sensor.  It  would  be  an  inordinate  task  to 
attempt  an  analysis  for  all  model  and  path  profiles. 
However,  the  results  are  based  on  the  absorption  from  CO^  in 
the  atmosphere,  and  the  concentration  of  CC^  is  relatively 
independent  of  input  conditions  (Ref  3:  V, 102-103). 

When  integrating  spectral  intensities  to  obtain  a  total 
band  intensity,  Simpson's  Rule  will  be  used  with  5  cm-* 
increments.  Although  the  resolution  of  current  spectro¬ 
meters  is  better  than  this,  the  LOWTRAN  transmittance  can 
only  be  calculated  in  5  cm-*  increments.  The  LOWTRAN 
transmittance  is  actually  an  average  over  a  20  cm-*  band. 
The  resolution  of  the  analysis  is  therefore  limited  by 
LOWTRAN. 


II.  Infrared  Signatures 

Source  Types 

There  are  numerous  ways  to  classify  sources  of  electro¬ 
magnetic  radiation.  Most  classification  schemes  rely  on  the 
spectral,  temporal,  or  spatial  properties  of  the  source  to 
determine  its  category.  In  order  to  predict  range  by  ana¬ 
lyzing  atmospheric  effects,  the  spectral  properties  of 
sources  are  the  most  important. 

In  general,  a  source  can  be  spectrally  categorized  as  a 
blackbody,  greybody,  or  selective  radiator.  A  blackbody  can 
be  defined  as  a  perfect  radiator,  or  one  that  radiates  the 
maximum  number  of  photons  that  any  body  in  thermodynamic 
equilibrium  at  the  same  temperature  can  radiate  (Ref  3:1,2). 
Other  sources  are  compared  to  blackbodies  by  a  quantity 
called  emissivity.  Emissivity  €  is  defined  as  the  ratio  of 
the  exitance  of  a  given  source  to  that  of  a  blackbody  (Ref 
3:1,28) : 


€(V) 


M(y,T) 

Mgg  (  J'jT) 


(1) 


where 

M(P,T)  -  spectral  exitance  of  a  source  at  temperature 
of  T 


M__(l>,T)  =  spectral  exitance  of  a  blackbody  at  same 
temperature 


A  source  whose  emissivity  is  a  constant  indepen¬ 
dent  of  wavelength  is  a  greybody,  while  sources  whose 


emissivity  is  a  function  of  wavelength  are  called  selective 
radiators.  In  many  instances,  selective  radiators  have 
emissivity  functions  that  cannot  be  easily  represented 
analytically.  Theoretical  model  calculations  can  produce 
tabulated  spectral  values,  but  the  spectral  signature  of 
these  sources  is  best  represented  by  measured  results. 

Blackbody  Sources.  When  a  cavity  is  an  isothermal 
enclosure  at  a  temperature  T,  a  small  hole  in  the  enclosure 
will  emit  blackbody  radiation  whose  spectral  exitance  can  be 
calculated  from  Planck's  Law  (Ref  2:126): 


MBB( 

where 


[W/cm^-cnfl  ] 


V  =*  wavenumber  in  cm”"1 
C±  *  2  7T he2  =  3.74  X  10‘12W-cm2 
C2  «  (hc/k)  =  1.439  cm-K 
Here,  c  is  the  velocity  of  light,  h  is  Planck's 


(2) 


constant, 


and  k  is  Boltzmann's  constant. 


The  spectra  for  several  blackbody  radiators  are  shown 
in  Figure  2.  If  the  temperature  of  a  blackbody  source  is 
known,  these  spectra  can  be  used  to  estimate  the  spectral 
band  over  which  a  sensor  should  operate  to  detect  the 
source. 

When  comparing  actual  sources  to  blackbody  sources, 
similar  units  must  be  used.  The  available  source  files  are 
recorded  in  units  of  radiant  intensity.  Assuming  a  small 


Figure  2.  Intensity  of  Blackbodies  at  Several 
Temperatures 

field-of-view  detector  and  normal  incidence,  a  point  source 
approximation  can  be  used  so  that  the  spectral  intensity  can 
be  calculated  from  the  spectral  exitance  by: 

I  (V)  «  (A /  7T  )  M  (V)  [W/(sr-cm_1)  ]  (3) 

e  e 

where 

1^ (i/)  =  spectral  radiant  intensity  at  the  source 
A  =  area  of  source 

To  remove  the  spectral  dependence,  the  spectral  intensity  is 
integrated  over  the  spectral  bandwidth  of  the  detector: 


I  =  f*2  I  ( v )  dv 

e  J  v*  e 


[W/sr] 


£  £ 
>>  ^ 


®! 


*  •  1  »  d 


where 

V  =  lower  wavenumber  of  bandwidth 

?2  =  upper  wavenumber  of  bandwidth 

Selective  Radiators.  Because  of  their  applicability  to 
infrared  ranging,  selective  radiators  resembling  aircraft 
jet  engines  will  be  used.  Aircraft  have  special  IR  signa¬ 
tures  resulting  from  the  volume  of  hot  gases  emitted  by  the 
engine  exhaust.  Missiles,  flares,  and  high-temperature 
industrial  smokestacks  also  exhaust  hot  gases  and  their 
spectra  therefore  resemble  that  of  aircraft.  An  example  of 
a  missile  spectrum  is  shown  in  Figure  3.  The  radiation 
resulting  from  the  vibrational  transitions  of  CO 2  and  H20 
are  apparent  between  2000-2500  cm"^  and  between  3000-4000 
cm”*.  If  the  spectral  resolution  of  the  detector  is  suffi¬ 
cient,  the  fine  structure  representing  transitions  between 
various  rotational  energy  states  of  the  molecules  composing 
the  exhaust  can  also  be  seen. 

For  hydrocarbon  fuels,  the  exhaust  consists  mainly  of 
H^O  vapor  and  CO^.  Because  of  these,  radiation  is  most 
intense  around  2300  cm-1  and  around  3700  cm-1.  However,  the 
atmosphere  is  also  composed  of  H^O  and  CO^  and  is  highly 
absorbent  in  these  spectral  bands.  Since  the  gases  are 
hotter  than  the  atmosphere,  gas  molecules  are  excited  to 
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Figure  3.  Typical  Missile  Spectrum  (ARI  #003) 


higher  rotational  energy  levels  than  the  atmospheric  mole¬ 
cules.  As  the  hot  gas  molecules  decay  to  the  ground  state, 
they  radiate  at  frequencies  that  are  outside  of  the  strong 
atmospheric  absorption,  and  the  radiation  propagates  with 
much  less  attenuation  (Ref  3:11,76-79).  This  effect  causes 
the  apparent  lack  of  radiation  in  the  band  between  2225  and 
2375  cm-*.  The  peak  of  radiation  remaining  near  2200  cm”* 
is  called  the  red  spike  while  the  peak  near  2400  cm”*  is 
called  the  blue  spike.  If  an  aircraft  jet  engine  is  running 
in  afterburner,  the  higher  temperature  causes  more  transi¬ 
tions  outside  of  the  atmospheric  absorption  region  to  occur, 
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Figure  4.  Aircraft  Spectrum  at  Different  Power  Settings 


and  the  red  spike  becomes  very  wide.  Because  the  blue  spike 
is  near  a  bandhead,  its  width  remains  relatively  constant. 
A  comparison  of  the  exhaust  spectra  for  aircraft  at  idle  and 
military  power  is  shown  in  Figure  4.  The  radiation  band  is 
obviously  wider  when  the  engine  is  in  military  power. 

In  addition  to  the  selective  radiation  caused  by  molec¬ 
ular  transitions,  the  IR  spectrum  contains  a  continuum  of 
radiation  caused  by  blackbody-type  sources  near  the  target. 
The  hot  tailpipe  and  outer  surfaces  of  an  aircraft,  or  the 
smokestack  of  a  refinery  flare,  all  contribute  to  the 
overall  IR  signature. 


Atmospherics 

The  spectral  signature  of  any  radiative  source  is 
altered  as  it  propagates  through  the  earth's  atmosphere. 
Radiation  is  attenuated  by  scattering  and  absorption 
processes  that  depend  upon  the  frequency  of  the  radiation. 
The  received  irradiance  is  proportional  to  the  initial 
radiance,  the  density  of  the  attenuating  medium,  and  the 
distance  traversed  (Ref  3:IV,3). 

Several  practical  methods  exist  for  calculating  the 
atmospheric  spectral  transmittance  to  accuracies  of  5  to  107. 
(Ref  3:V,3).  The  LOWTRAN  method  uses  a  functional  form  of 
the  transmittance  that  is  determined  by  a  fit  of  a  special¬ 
ized  curve  to  a  plot  of  actual  data  (Ref  3:V,47).  The 
results  have  been  computerized  and  the  FORTRAN  code  is 
available  from  the  Air  Force  Geophysics  Laboratory  (Ref  1). 
LOWTRAN  is  frequently  updated  and  the  version  used  in  this 
analysis  is  called  LOWTRAN  5. 

The  LOWTRAN  code  calculates  transmittance  averaged  over 
20  cm-*  intervals  in  steps  of  5  cm-*  from  350  to  40,000  cm-* 
(Ref  1:14).  The  1962  U.S.  Standard  atmosphere  and  five 
seasonal  models  are  provided  as  input  options  for  running 
the  code.  In  addition,  user-defined  or  measured  data  may  be 
inserted  to  represent  any  atmospheric  conditions. 

To  represent  the  irradiance  detected  by  an  infrared 
sensor,  computer  data  files  containing  the  transmittance  for 
an  atmospheric  model  at  several  ranges  were  created.  The 
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files  are  generated  only  for  the  frequencies  over  which  the 
analysis  is  anticipated.  An  array  that  contains  the  inten¬ 
sity  of  a  source  can  then  be  multiplied  by  the  appropriate 
transmittance  in  the  data  file  and  divided  by  the  square  of 
the  range  to  obtain  irradiance  at  the  sensor. 


Background  Contributions 

In  the  infrared  region  of  the  electromagnetic  spectrum, 
background  radiation  has  a  significant  effect  on  the  amount 
of  radiation  seen  by  a  sensor.  Most  detection  systems  are 
contrast  sensors  and  the  primary  response  generated  from  a 
source  is  due  to  the  contrast  within  the  scene.  Any  black- 
body  radiation  from  the  background  is  also  sensed  and  be¬ 
comes  a  part  of  the  total  measured  spectrum.  If  a  target- 
plus-background  spectrum  is  obtained,  the  target-only  spec¬ 
trum  could  be  calculated  by  spectrally  subtracting  the 
background-only  spectrum.  When  using  this  technique,  care 
must  be  exercised  so  that  the  background  spectrum  is 
obtained  under  the  same  conditions  that  the  composite  spec¬ 
trum  was  obtained  (Ref  4:1,1). 

In  addition  to  the  steady  background  radiation,  a  clut¬ 
tered  background  such  as  clouds  or  terrain  that  scatters 
radiation  exhibits  an  intensity  pattern  that  varies  from 
point  to  point.  A  sensor  scanning  such  a  background  gener¬ 
ates  a  response  that  appears  as  noise  on  the  signal.  In 
many  applications,  this  clutter  noise  limits  the  detection 


Computer  Files 

Two  sets  of  data  files  designed  to  be  used  for  spectral 
discrimination  analyses  were  available  from  previous  con¬ 
tracts.  Aerodyne  Research  Inc.  has  supplied  spectra  that 
are  computer  generated,  and  validated  through  comparison 
with  actual  measurements.  The  Air  Force  Geophysics  Labora¬ 
tory  has  provided  spectra  that  are  taken  from  actual  field 
measurements. 

ARI  Files .  The  library  of  spectra  obtained  from 
Aerodyne  are  referred  to  as  ARI  files.  Although  numerous 
source  spectra  were  available,  only  nine  represented  targets 
that  resembled  aircraft  and  missiles.  One  of  these  files  is 
measured  data  from  a  missile  firing,  and  the  remainder  are 
computer  simulations  of  flares  typically  seen  near  refin¬ 
eries.  Different  gas  and  particle  concentrations  cause 
minor  variations  in  the  spectra,  but  all  are  similar. 

Table  I  is  a  list  of  the  files  that  were  used  with  a 
short  description  of  each.  Several  of  the  spectra  are 
plotted  in  Figure  5.  Computer  generated  files  are  charac¬ 
terized  by  no  noise  noticeable  in  the  areas  of  low  intensity 
radiation.  Except  for  the  missile  spectrum,  all  of  the  ARI 
files  have  no  noise  associated  with  the  signals.  The  mis¬ 
sile  file  is  actually  measured  data  and  not  computer  gener¬ 
ated,  so  it  does  have  noise  associated  with  the  signal. 

The  range  of  ARI  files  is  simulated  by  multiplying  each 
wavenumber  by  a  corresponding  LGWTRAN  transmittance.  Using 


File  name 

Temp 

Description 

ARILIB . 003 

Unknown 

Missile  firing 

ARILIB.015 

1800  K 

Gas  burnoff  flare 

ARILIB. 017 

1800  K 

Gas  burnoff  flare 

ARILIB. 019 

1800  K 

Gas  burnoff  flare 

ARILIB. 035 

1600  K 

Arc  flare 

ARILIB. 050 

1800  K 

Gas  burnoff  flare 

ARILIB. 052 

1800  K 

Gas  burnoff  flare 

ARILIB. 053 

1800  K 

Gas  burnoff  flare 

ARILIB. 054 

1800  K 

Gas  burnoff  flare 
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Figure  5b.  ARI  Flare  Spectrum  (#017) 


this  method,  sources  representing  several  IR  targets  can  be 
simulated  with  different  atmospheric  conditions  and  at  dif¬ 
ferent  ranges. 

AFG  Files .  Actual  measured  spectra  obtained  from  the 
Air  Force  Geophysics  Laboratory  are  referred  to  as  AFG 
files.  These  measurements  were  obtained  from  an  aircraft 
flying  over  industrial  areas  of  the  South  and  Midwest  in 
June  1978  with  a  one  degree  f ield-of-view  spectrometer  (Ref 
5:100).  The  field-of-view  projection  is  larger  than  the 
diameter  of  most  of  the  sources,  so  the  original  spectra 
contain  contributions  from  the  backgrounds.  Background 
spectra  near  the  targets  were  also  recorded,  and  were 
subtracted  from  the  target  spectra  to  provide  source-only 
spectra.  Because  the  ranges  and  background  radiances  are 
somewhat  uncertain,  the  spectra  have  associated 
uncertainties  of  perhaps  +257»  (Ref  5:106).  Unfortunately, 
the  ranges  of  all  sources  are  only  between  2.6  and  4.7  km. 

Several  of  the  sources  do  not  represent  hot  gas 
sources,  but  are  more  characteristic  of  cool  gases  and 
greybody  emitters.  Cool  gases  are  sources  that  have  only  a 
small  amount  of  radiation  that  is  emitted  beyond  the 
atmospheric  absorption  band  (<1000  K);  hot  gases  emit  well 
beyond  the  atmospheric  absorption  band.  The  actual  type  of 
these  sources  falls  into  one  of  the  following  categories: 

1)  Coke  oven/blast  furnace 

2)  Smoke  and  flare  stacks 


Table  II 


Air  Force  Geophysics  Lab  Files 


File  name 

Type 

Temp 

Location 

Range 

AFGLIB.001 

Hot  Ore 

819  K 

Indiana 

3.5 

AFGLIB. 002 

Hot  Slag 

845  K 

Indiana 

3.5 

AFGLIB. 003 

Coke  Oven 

496  K 

Indiana 

3.9 

AFGLIB. 004 

Flare  Stack 

None 

Indiana 

4.7 

AFGLIB. 005 

Mill  Stacks 

600  K 

Indiana 

4.4 

AFGLIB. 006 

Blast  Furn. 

600  K 

Michigan 

4.7 

AFGLIB. 007 

Coke  Oven 

481  K 

Michigan 

2.8 

AFGLIB. 008 

Stacks 

800  K 

Michigan 

4.4 

AFGLIB. 009 

Hot  Metal 

800  K 

Michigan 

3.2 

AFGLIB. 010 

Stacks 

600  K 

Michigan 

3.2 

AFGLIB. Oil 

Flare  Stack 

1700  K 

Kentucky 

3.3 

AFGLIB. 012 

Hot  Bldgs. 

500  K 

Kentucky 

2.9 

AFGLIB. 013 

Stacks 

400  K 

Alabama 

2.9 

AFGLIB. 014 

Hot  Bldgs. 

600  K 

Alabama 

3.1 

AFGLIB. 015 

Hot  Metal 

1500  K 

Alabama 

2.6 

AFGLIB. 016 

Hot  Bldgs 

575  K 

Kentucky 

2.7 

3)  Hot  metals  and  ore 

4)  Hot  buildings 

One  of  each  of  these  types  of  source  spectra  is  plotted  in 
Figure  6,  and  a  list  of  all  the  files  is  given  in  Table  II. 
The  temperature  listed  in  Table  II  is  the  temperature  of  a 
blackbody  source  whose  spectrum  closest  resembles  that  of 
the  actual  source.  It  is  not  necessarily  the  actual  source 
temperature  because  the  source's  spatial  distribution  is  not 
considered.  Because  the  spectrometer  had  a  high  resolu¬ 
tion,  the  fine  structure  of  the  emission  is  noticeable  in 
Figure  6.  This  will  later  cause  problems  because  the  ARI 
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Fig  6c.  AFG  Smoke  Stack  Spectrum  (#010) 


files  multiplied  by  LOWTRAN  generate  low  resolution  spectra 
in  which  none  of  the  atmospheric  fine  structure  is  present. 

Because  of  their  security  classification,  actual 
sources  representing  aircraft  and  missile  signatures  could 
not  be  obtained.  However,  the  AFG  files  do  represent  actual 
measured  results  and  are  therefore  used  to  test  the  ranging 
algorithm.  Unfortunately,  the  range  of  these  sources  only 
varies  from  2.6  to  4.7  km,  and  thus  they  provide  only  a 
limited  check  of  the  algorithm. 


m  m 
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III.  Spectral  Discrimination 


Pattern  Recognition 

Spectral  pattern  recognition  is  the  process  of  identi¬ 
fying  a  class  of  events  by  analyzing  the  spectrum  that  has 
been  produced  or  influenced  by  the  event  (Ref  6:2699). 
Ideally  the  process  will  correctly  identify  all  of  the 
desired  events  and  reject  all  background,  interference,  and 
noise. 

The  received  spectrum  of  any  source  may  be  represented 
by  N  spectral  amplitudes  (S ,s2, • • • ,SN)  at  the  corresponding 
N  different  wavelengths  or  colors.  The  spectrum  can  then  be 
represented  in  an  N-dimensional  color  space  by  the  coordi¬ 
nates  (S^ jS^, . . . ,S^)  along  N  orthogonal  axes.  A  set  of 
spectra  from  one  class  of  events  can  be  represented  by  a 
group  of  points  in  the  color  space.  Assuming  that  the 
spectra  from  the  class  are  similar,  the  points  will  form  a 
cluster. 

If  a  different  class  is  represented  by  the  same  scheme, 
two  clusters  will  form  in  color  space.  If  little  or  no 
overlap  exists  between  the  two  clusters,  the  classes  can  be 
accurately  distinguished.  Figure  7  shows  an  example  of  two 
classes  of  events  plotted  in  a  two-dimensional  color  space. 

Because  the  spectrum  of  any  event  or  class  of  events  is 
not  unique,  pattern  recognition  becomes  a  statistical  prob¬ 
lem.  To  simplify  the  statistical  calculations,  normal 
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Figure  7.  Spectral  Classes  in  Color  Space  (Ref  6) 


Gaussian  statistics  are  assumed.  Although  the  data  obtained 
from  infrared  spectra  does  not  entirely  support  normal 
statistics,  it  has  been  shown  that  a  decision  rule  based  on 
normal  distributions  performs  nearly  as  well  as  rules 
derived  from  the  empirical  frequency  of  the  data  (Ref 
7:162).  By  using  normal  statistics,  calculations  are 
simplified  and  the  distributions  can  be  uniquely  character¬ 
ized  by  their  mean  vector  and  covariance  matrix  (Ref  8:22). 

Distributions  in  color  space  can  be  partitioned  into 
separate  regions  called  decision  regions.  One  decision 
region  must  be  associated  with  each  class  of  events  or 
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Figure  8.  Decision  Boundaries  for  Normal  Distributions 
(Ref  8:53) 


cluster,  and  should  contain  most  of  the  points  for  that 
cluster.  Any  new  point  from  an  unclassified  source  can  be 
assigned  according  to  the  decision  region  into  which  it 
falls. 

There  several  methods  available  for  determining  the 
decision  boundaries.  Using  a  maximum  likelihood  criteria 
(Ref  8:53),  construction  of  the  decision  regions  is  easily 
visualized  by  plotting  contours  of  equal  density  as  in 
Figure  8.  The  decision  function  that  separates  the  two 
distributions  is  defined  by  the  intersections  of  contours  of 
equal  density.  For  Gaussian  distributions,  the  contours  are 
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ellipses  and  the  desision  functions  are  quadratic.  When  the 
two  distributions  have  equal  covariance,  the  decision 
functions  become  linear.  A  method  to  calculate  the  slope  of 
the  line  separating  normal  distributions  of  equal  covariance 
is  discussed  in  Appendix  A. 

Gaussian  Statistics 

A  vector  whose  components  are  random  variables  is 
called  a  random  vector.  A  random  vector  is  fully  character¬ 
ized  by  its  distribution  or  density  function.  Often  these 
functions  cannot  be  easily  determined,  or  are  too  complex  to 
be  of  practical  use.  It  is  therefore  convenient  to  use  less 
accurate  but  more  computable  density  functions  to  represent 
some  types  of  distributions. 

When  using  Gaussian  statistics,  the  density  function  of 
a  multivariate  distribution  is  completely  represented  by  its 
mean  or  expected  vector  and  its  covariance  matrix.  The  mean 
vector  of  a  random  vector  X  is  defined  as  (Ref  8:16): 

H  =  E[X]  = 

where 

E[X]  =  expectation  value  of  vector  X 
p(X)  =  density  function  of  X 
s  =  all  X  space 

The  covariance  matrix  C  is  defined  by  (Ref  8:17; 10 : 34 ) : 


L 


X  p(X)  dX 


(5) 


C  =  E[ (X-M)  ( X-M) ] 


=  (X-H)T(X-M)  p(X)  dX  (6b) 
Js 

where  T  represents  the  transpose  of  a  vector.  The  diagonal 
components  of  this  matrix  are  the  variances  of  individual 
random  variables  and  the  off-diagonal  components  are  the 
covariances  of  the  two  corresponding  variables. 

When  only  a  finite  number  n  of  spectra  are  available  to 
characterize  a  distribution,  the  mean  vector  is  calculated 
by  (Ref  9:236) 


M  =  (1/n)  y*. X 


The  covariance  matrix  is  estimated  by 


C  =  [l/(n-l)  ]  ^(X-fi)T(X-M) 


Sometimes  it  is  convenient  to  normalize  covariance 
matrices  by  converting  the  individual  covariance  terms  to 
correlation  coefficients  r^  as 


=  *ii  *jj 


where 

c^j  =  matrix  element  of  C 
aii  =  stant*arc*  deviation  of  band  i 
Because  the  variance  is  the  square  of  the  standard 


deviation,  the  diagonal  elements  of  the  correlation  matrix  R 
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are  unity. 

When  using  Gaussian  statistics,  the  distribution  func¬ 
tion  pOO  is  given  in  terms  of  the  mean  vector  and  covar¬ 
iance  matrix  by  (Ref  10:35): 

pCX)  =  (2  n  }“d/2  |  C  | — 1  ^ 2  exp [ h. ( X-H ) TC-1  ( X-M )  ]  (10) 

where 

d  =  dimension  of  X 
|C|  =  determinant  of  C 
C“*  =  inverse  of  C 

Cluster  Discrimination 

To  reduce  the  probability  of  an  erroneous 
classification,  the  separation  between  distributions  must  be 
large.  The  Bhattacharyya  distance  is  one  measure  of  this 
separation  and  is  defined  by  (Ref  8:268): 

B  =  ”ln  Jj^Pl(X)  P2(*}^  dX  (11) 

By  assuming  normal  statistics,  the  Bhattacharyya  distance 
can  be  written  as  (Ref  8:268): 

B  =  (1/16) (H1-H2)T(C1+C2)"1 (Hj-Hj) 

+  ^ln[^(C1+C2)  ICir*IC2r*]  (12) 

This  distance  is  zero  when  two  distributions  are  identical, 
and  is  infinite  when  there  is  no  overlap. 

The  Bhattacharyya  distance  is  the  sum  of  two  terms. 
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The  first  term  in  the  equation  is  called  the  Mahalanobis 
distance  and  is  essentially  the  square  of  the  distance 
between  the  means  of  each  distribution.  The  distance  is 
measured  relative  to  the  average  of  the  covariances.  The 
second  term  is  an  added  separation  arising  from  different 
covariances.  If  the  means  of  the  distributions  are  equal, 
the  Mahalanobis  distance  vanishes,  but  there  may  still  be 
separation  from  the  second  term.  If  the  covariances  are 
equal,  the  second  term  vanishes  and  the  only  separation  is 
from  the  difference  of  the  means. 

Using  a  maximum  likelihood  criterion,  the  probability 
of  an  erroneous  classification  can  be  approximated  in  terms 
of  the  Bhattacharyya  by  the  expression  (Ref  11) 

E  =  ( 2B2 JT  )"^exp(-B2)  (13) 

Error  for  a  Bhattacharyya  distance  of  1.0  is  calculated  to 
be  approximately  0.147  while  a  Bhattacharyya  distance  of  2.0 
results  in  an  error  of  only  0.004.  For  Gaussian  statistics, 
a  Bhattacharyya  separation  of  2.0  therefore  represents  small 
overlap  between  two  distributions. 

In  some  applications  of  pattern  recognition,  it  is 
practical  to  use  only  a  few  wavelengths  or  spectral  bands 
when  analysing  source  spectra.  Optimum  bands  are  found  by 
changing  the  combination  of  wavelengths  until  the  Bhatta¬ 
charyya  distance  between  distributions  in  color  space  is 
maximized.  The  least  number  of  bands  is  sought  that  still 


yvvv 


permits  adequate  separation  between  distributions. 
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Intensity  Normalization 

The  intensity  of  any  source  has  a  direct  effect  on  the 
location  of  the  point  plotted  in  color  space.  Bright 
sources  are  located  farther  from  the  origin  than  dim  sources 
of  the  same  class  because  the  spectral  amplitude  of  all 
bands  will  be  greater  for  the  bright  source.  When  trying  to 
discriminate  certain  classes  of  events  by  spectral  pattern 
recognition,  it  is  sometimes  desirable  to  obtain  results 
that  are  independent  of  source  intensity. 

In  the  development  of  spectral  representation  using 
color  space,  a  point  representing  one  spectrum  of  a  class  of 
events  was  plotted  using  rectangular  coordinates 
(S^ ,S2 , . . . ,SN)  defined  by  the  spectral  amplitude  in  N  bands. 


:  S 

Using  only  three  bands  for  comparison,  any 

point 

in  color 

I  c 

"  - 

^  ", 

space  can  also  be  represented  with  spherical 

coordinates 

(p,e,4>)  by  (Ref  12:385): 

p  —  ( +^2  ' 

(14a) 

4>  » 

e  =  arctan(S2/S^ ) 

(14b) 

2  2 

4>  =  arccos[ (S^) / (S^  +S^ 

+S^] 

(14c) 

Using  spherical  coordinates, 

the  intensity 

of  a 

source  is 

indicated  by  the  magnitude  of  p,  and  the  relative  inten¬ 
sities  among  the  three  bands  are  given  by  the  values  of  e 
and  <j>.  If  it  is  desirable  to  have  a  color  space  represen¬ 
tation  that  is  independent  of  source  intensity,  the  point 
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can  be  plotted  using  a  two-dimensional  coordinate  system 
whose  axes  are  represented  in  rectangular  coordinates  by  the 
values  of  e  and  4>. 

If  only  two  bands  are  used  in  the  color  space  represen¬ 
tation,  the  point  can  similarly  be  located  using  plane  polar 
coordinates  where  (Ref  12:378): 


r  *  (S?  +  Sl)h  (15a) 
e  =  arctan(S2/S1)  (15b) 

If  a  color  space  representation  independent  of  source  inten¬ 
sity  is  required,  the  two-dimensional  point  can  be  repre¬ 
sented  in  only  one  dimension  defined  by  the  value  of  e.  In 
addition  to  removing  the  dependence  of  the  distributions  on 
source  intensities,  this  "reduced"  color  space  represen¬ 
tation  simplifies  the  calculation  of  the  statistical  mea¬ 
sures  by  reducing  the  vectors  and  matrices  by  one  dimension. 


IV.  Spectral  Optimization 


Band  Selection 

It  has  already  been  mentioned  that  blackbody  curves  are 
useful  to  predict  the  spectral  band  over  which  a  sensor  must 
operate.  In  Figure  2,  blackbody  curves  were  plotted  for 
several  realistic  temperatures.  Radiation  for  blackbody 
sources  near  1000  K  is  most  intense  between  1000  and  4000 
cm”1. 

Actual  sources,  however,  do  not  emit  broadband  radia¬ 
tion  as  blackbodies  do.  An  example  of  a  missile  spectrum 
was  plotted  in  Figure  3.  Sensors  that  are  to  detect  this 
type  of  source  can  be  further  limited  to  the  regions  between 
2000  and  2500  cm-1  or  3300  and  3800  cm”1.  Much  of  this 
radiation  is  absorbed  by  the  atmosphere,  though.  Because 
the  transmittance  of  the  band  between  3300  and  3800  cm”1  is 
lower,  and  scattered  solar  radiation  is  more  pronounced,  the 
optimum  spectral  region  is  between  2000  and  2500  cm”1. 

To  avoid  range  errors  that  may  be  caused  by  climatic 
differences,  spectral  bands  should  be  chosen  that  minimize 
the  variation  in  transmittance  between  different  atmospheric 
models.  L0WTRAN  transmittance  for  several  models  at  the 
same  range  is  plotted  in  Figure  9.  Variations  from  H^O 
vapor  become  very  pronounced  below  2175  cm”1.  Wavenumbers 
below  this  value  will  not  be  used.  Transmission  in  the  CO^ 
band  between  2275  and  2375  cm”1  is  negligible  for  all 


Figure  9.  LOWTRAN  Model  Comparison  at  2  km 

models,  so  this  band  cannot  be  used  either.  Above 
2375  cm--*-,  a  sensor  would  be  detecting  a  source's  blue 
spike.  Because  the  blue  spike  is  very  narrow  and  source 
spectra  in  this  region  vary  considerably,  the  wavenumbers 
above  2375  cm--*-  will  also  be  avoided.  The  only  band  that 
remains  is  between  2175  cm-^  and  2275  cm-^.  The  atmospheric 
transmittance  at  several  ranges  is  plotted  in  Figure  10  for 
this  spectral  region.  The  regular  increase  in  CO2  absorp¬ 
tion  is  the  atmospheric  effect  that  will  be  used  to  deter¬ 


mine  the  range  of  a  source  by  examining  its  spectrum. 

One  of  the  limitations  of  LOWTRAN  is  the  spectral 
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Figure  10.  LOWTRAN  Range  Comparison  for  Model  6 

resolution  of  20  cm-1-.  Minimum  bandwidths  will  be  restric¬ 
ted  to  25  cm-1  because  of  this  limitation.  If  smaller  band- 
widths  are  required  to  resolve  spectral  differences,  a  more 
accurate  atmospheric  model  will  have  to  be  used  for  the 
analysis . 

Bhattacharyya  Measure 

The  spectral  region  between  2175  and  2275  cm-*  was 
divided  into  four  bands  as  follows: 


band 

1: 

2175-2200  cm-1 

band 

2: 

2200-2225 

cm-* 

band 

3: 

2225-2250  cm-1 

band 

4: 

2250-2275 

cm-* 
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By  choosing  any  combination  of  these  bands,  the  spectrum  of 
any  source  can  be  represented  in  color  space.  Initially 
only  two  bands  were  used  for  comparison,  so  color  space  is  a 
two-dimensional  coordinate  system  whose  axes  are  the  inten¬ 
sities  in  each  band.  To  remove  the  dependence  of  the  color 
space  representation  on  the  intensity  of  the  source,  the 
spectrum  was  represented  in  reduced  color  space.  For  two 
bands,  reduced  color  space  is  one  dimension  whose  value  is 
the  arctangent  of  the  ratio  of  the  intensity  in  the  two 
bands . 

The  first  sources  used  were  blackbody  spectra  multi¬ 
plied  by  LOWTRAN  transmittance.  Several  blackbody  functions 
at  temperatures  between  500  K  and  2000  K  were  multiplied  by 
three  different  LOWTRAN  models.  The  models  chosen  represent 
midiattitude  summer  (model  2),  midiattitude  winter  (model 
3),  and  the  1962  U.S.  Standard  (model  6).  Ranges  between  2 
km  and  16  km  were  calculated.  The  resulting  spectra  were 
then  integrated  over  the  selected  frequency  bands  to  deter¬ 
mine  the  location  of  the  spectral  representation  in  reduced 
color  space. 

All  the  source  combinations  for  one  range  form  a 
distribution  of  points  in  reduced  color  space.  The  mean  and 
covariance  matrix  for  the  distribution  were  calculated.  In 
one  dimension,  the  covariance  matrix  is  simply  the  square  of 
the  standard  deviation  of  the  angle  e.  The  distributions 
for  all  the  ranges  used  were  calculated  using  this  method. 
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To  determine  the  separation  between  each  range  distri¬ 
bution,  the  Bhattacharyya  distance  was  calculated.  Large 
values  for  the  Bhattacharyya  distance  mean  that  the  range 
distributions  are  well  separated.  Obviously,  the  actual 
distance  will  depend  on  the  range  increment  used.  Because 
the  Bhattacharyya  distance  is  only  used  as  a  comparison 
between  band  choices,  the  range  increment  is  not  critical. 
It  does  however  give  an  idea  of  the  resolution  that  may  be 
obtained  using  this  method.  If  the  Bhattacharyya  separation 
between  2  km  range  distributions  is  large,  then  2  km  range 
resolution  may  be  possible.  If  there  is  large  overlap 
between  2  km  distributions,  there  will  be  large  uncertainty 
as  to  which  distribution  an  unclassified  source  should  be 
assigned. 

After  calculating  the  Bhattacharyya  separation  between 
each  range  distribution,  the  band  choice  was  changed  and  new 
distributions  were  formed.  Bhattacharyya  distances  between 
the  new  distributions  were  calculated  and  compared  with  the 
results  of  the  previous  bands.  The  band  choice  that  pro¬ 
vides  the  highest  separation  between  range  distributions 
will  provide  the  best  resolution  when  trying  to  assign  un¬ 
classified  sources  to  a  specific  range. 

In  Figure  11,  the  Bhattacharyya  distances  obtained  with 
several  possible  band  combinations  are  plotted.  Band  4 
combinations  are  not  shown  because  the  Bhattacharyya  dis¬ 
tances  obtained  were  small.  The  intensity  in  this  band  was 
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Figure  11.  Bhattacharyya  Separation  for  Blackbody 
Sources  Using  Two  Colors 


too  low  at  all  ranges  to  allow  sufficient  separation  between 
range  clusters.  At  close  ranges,  the  separation  is  very 
good  for  all  other  band  choices.  As  range  increases 
however,  the  separation  drops  rapidly.  This  is  because  the 
intensity  of  the  highest  wavenumber  band  becomes  very  small 
due  to  C(>2  absorption,  and  the  ratio  of  the  two  bands  ap¬ 
proaches  0.  As  the  ratio  approaches  0,  so  does  the  angle  e 
that  defines  the  location  of  the  distributions  in  reduced 
color  space.  The  range  distributions  tend  to  congregate 
near  e  «  0,  and  the  Bhattacharyya  distances  therefore  become 


The  optimum  bands  for  blackbody  sources  is  band  1 
(2175-2200  cm”*)  and  band  2  (2200-2225  cm”*).  Bhattacharyya 
distances  between  range  clusters  remains  above  2.0  until 
10  km  using  only  two  colors  for  comparison.  Although  actual 
sources  are  not  well  represented  by  blackbody  functions, 
this  analysis  has  shown  that  IR  signatures  can  be  used  to 
classify  sources  according  to  range. 

Source  files  in  the  ARI  library  were  analysed  in  the 
same  manner.  Because  of  the  number  of  source  files  avail¬ 
able  and  the  computer  time  required  to  read  and  reduce  the 
data  files,  LOWTRAN  was  limited  to  the  1962  U.S.  Standard 
model  atmosphere  and  23  km  haze  visibility.  Deviations 
caused  by  atmospheric  variations  are  discussed  in  the 
results.  Again,  only  two  bands  were  used  to  begin  the 
analysis. 

The  Bhattacharyya  distances  are  plotted  in  Figure  12 
for  all  possible  band  combinations.  The  optimum  separation 
is  obtained  with  band  2  (2200-2225  cm”*)  and  band  3  ( 2225— 
2250  cm”*).  With  these  bands,  the  Bhattacharyya  distance 
remains  above  2.0  only  through  5  km.  This  separation  is  not 
sufficient  to  attempt  a  ranging  algorithm,  so  three  bands 
for  comparison  will  have  to  be  used. 

It  should  be  noted  that  these  bands  are  at  a  higher 
wavenumber  than  the  bands  obtained  with  blackbodies.  With 
blackbodies,  it  was  found  that  the  absorption  near  the  CO^ 
band  became  too  large  for  higher  wavenumber  bands.  Best 


Figure  12.  Bhattacharyya  Separation  for  ARI  Sources 
Using  Two  Colors 

results  were  obtained  by  choosing  bands  as  far  away  from  the 
C(>2  band  as  possible.  With  actual  sources  however,  the 
radiation  intensity  begins  to  decrease  as  the  wavenumber 
moves  away  from  the  CC>2  band.  This  is  typical  of  cool  gas 
sources,  so  frequency  bands  must  be  chosen  adjacent  to  the 
C02  band  for  sources  that  are  cool  gas  radiators. 

Several  ARI  files  that  are  not  listed  could  not  be  used 
because  of  this  problem.  The  sources  represented  power 
plant  buildings  at  a  temperature  of  400  K,  and  the  spectra 
resembled  cool  gases.  When  the  spectral  representations  of 
these  sources  were  plotted  in  reduced  color  space,  the 


points  did  not  fall  near  the  proper  range  distribution. 
Cool  gas  sources  will  therefore  not  provide  accurate  range 
estimates  with  the  algorithm  that  will  be  determined. 

To  improve  the  Bhattacharyya  separation  between  range 
distributions,  the  ARI  spectra  were  analyzed  using  three 
bands.  The  only  choice  available  for  the  third  band  was 
between  band  1  and  band  4.  Because  the  intensity  of  band  4 
is  extremely  low  at  longer  ranges,  band  1  was  chosen  as  the 
third  band. 

The  integrated  intensity  in  all  three  bands  was  cal¬ 
culated  for  each  ARI  source  file  at  ranges  between  2  and  18 
km.  Using  three  colors,  the  spectral  representation  in 
reduced  color  space  is  a  point  in  a  two  dimensional  coordi¬ 
nate  system  whose  axes  are  the  values  of  the  spherical  coor¬ 
dinate  angles  e  and  <J>.  In  Appendix  B,  the  value  of  each 
band  intensity,  and  the  values  of  e  and  <J>,  are  tabulated  for 
all  sources  at  all  ranges.  In  addition,  the  mean  and  stan¬ 
dard  deviation  of  e  and  <J>  are  given  to  define  the  Gaussian 
distribution  parameters  for  each  range  cluster. 

In  Figure  13,  the  spectral  representations  for  all 
sources  are  plotted  in  reduced  color  space,  and  ellipses  are 
drawn  around  points  representing  sources  at  the  same  range. 
By  looking  at  the  graph,  it  is  obvious  that  there  is  good 
separation  between  all  of  the  range  clusters,  even  at  longer 
ranges.  The  Bhattacharrya  distances  between  the  distribu¬ 
tions  are  plotted  in  Figure  14.  Separation  increases  with 


range  until  13  km,  and  then  drops  rapidly.  Even  though  the 
distributions  appear  to  get  closer  together,  Bhattacharrya 
distance  is  measured  relative  to  the  variance  of  the  distri¬ 
butions  which  also  decreases.  The  net  effect  is  an  increase 
in  the  Bhattacharrya  distance  until  13  km  where  the  separ¬ 
ation  between  the  means  becomes  small  compared  to  the 
variance. 

Significant  improvement  is  obtained  in  Bhattacharyya 
separation  between  distributions  by  using  the  third  band. 
In  the  next  chapter,  the  mean  and  covariance  matrix  of  these 
distributions  will  be  used  to  form  the  decision  regions  into 
which  unclassified  sources  will  be  assigned.  Before  pro¬ 
ceeding,  the  effect  of  background  radiation  on  the  Bhatta¬ 
charyya  separation  will  be  investigated. 

Background  Fluctuations 

The  normalized  distributions  that  represent  range 
clusters  in  color  space  were  formed  assuming  no  background 
radiation  contributions.  Because  sensors  can  be  AC-coupled 
to  remove  the  DC  component  of  the  noise  signal,  the  only 
effect  a  background  will  have  is  to  add  fluctuations  to  the 
received  source  signals.  The  location  of  the  means  of  the 
range  clusters  in  color  space  will  therefore  remain  the 
same,  but  the  variances  and  covariances  will  increase. 
Because  of  this  effect,  the  Bhattacharyya  distances  previ¬ 
ously  obtained  represent  the  upper  limit  of  a  noise-free 
background.  Bhattacharyya  distances  will  decrease  with  the 


addition  of  noise  fluctuations. 

To  quantitatively  account  for  noise,  a  covariance 
matrix  representing  the  fluctuations  of  the  received  signals 
must  be  determined.  Because  the  range  covariance  is  statis¬ 
tically  independent  of  the  noise  covariance,  the  total 
covariance  matrix  is  the  sum  of  the  two  independent  covar¬ 
iance  matrices  (Ref  13:37). 

The  Environmental  Research  Institute  of  Michigan  (ERIM) 
has  published  comprehensive  and  well  documented  background 
data  in  the  mid-infrared  region  of  the  spectrum  (Ref  14), 
and  the  measured  data  has  been  extensively  analysed  in  an 
unpublished  report  (Ref  11).  In  Appendix  C,  an  expression 
for  the  noise  variance  in  the  bands  previously  used  is 
obtained  based  on  this  report.  The  value  obtained  for  the 
noise  variance  in  each  25  cm-1  band  is 

o2  =  3.75  x  10"3  R2  (W/sr)2  (16) 
where 

R  =  distance  from  source  to  sensor 

Because  the  range  clusters  are  formed  in  reduced  color 
space,  the  variance  in  each  band  needs  to  be  converted  to  a 
variance  in  terms  of  e  and  <f>.  This  is  done  through  differ¬ 
entiation  of  the  coordinate  transforms  and  the  details  are 
shown  in  Appendix  D.  The  result  is  an  expression  for  the 
variance  of  e  and  <j>  in  terms  of  the  variance  obtained  in 
Appendix  C.  To  complete  the  covariance  matrix,  a  value  for 


the  correlation  must  be  obtained.  Relying  on  the  measured 
data  by  ERIM,  the  correlation  is  estimated  to  be  0.7  (Ref 
11).  An  accurate  estimate  of  the  correlatioi.  s  not  pos¬ 
sible,  but  the  effect  on  the  results  is  minimal  anyway.  The 
diagonal  elements  of  the  covariance  matrix  (the  variances) 
have  the  most  significant  effect.  The  background  covariance 
matrix  that  was  added  to  the  range  distribution  covariance 
matrices  is  therefore: 


Cn  " 


(0.7)  (de)  (d<j>) 


(0.7)  (de)  (dd>) 


where  de  and  d<j>  are  defined  in  Appendix  D  by  equations  36b 
and  37c. 

The  covariance  matrix  for  noise  fluctuations  depends  on 
range  (de  and  d<j>)  because  the  area  of  the  background  seen  by 
the  sensor  depends  on  how  far  away  the  sensor  is  from  the 
background.  As  the  detected  area  increases,  so  does  the 
magnitude  of  the  fluctuations  in  the  received  signal.  Also 
associated  with  the  range  is  the  attenuation  of  the  fluctu¬ 
ations  because  of  atmospheric  absorption. 

The  covariance  matrices  of  the  range  distributions  were 
added  to  the  covariance  matrices  of  the  noise  fluctuations 
and  the  Bhattacharyya  distances  were  recalculated.  The 
effect  of  the  noise  covariance  on  the  Bhattacharyya  dis¬ 
tances  depends  on  the  source  intensity.  In  rectangular 
coordinates,  the  magnitude  of  the  fluctuations  is  constant, 
but  by  using  spherical  coordinates  to  locate  the  distri- 


Figure  15.  Bhattacharyya  Separation  for  ARI  Sources 
With  Noise  Added  (I  in  W/sr  at  2  km  in 
band  2200-2225)  a 

butions,  the  dependence  on  the  source  intensity  is  intro¬ 
duced.  Weak  sources  relative  to  the  background  signal  are 
affected  more,  so  the  separation  between  range  distributions 
decreases  more. 

The  results  of  the  calculations  are  plotted  in  Figure 
15.  The  plot  labeled  as  infinite  intensity  is  actually  the 
original  Bhattacharyya  measures  with  no  noise  added  (the 
result  is  the  same).  The  intensity  I  is  an  apparent 
intensity  at  2  km  in  the  band  2200-2225  cm-^.  The  apparent 
intensity  at  2  km  is  given  by: 
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where 


I  =  intensity  at  the  source 
t  =  transmittan-e  of  the  atmosphere 
I  is  called  an  apparent  intensity  because  it  does  not 
include  the  inverse  square  variation  with  range. 

As  expected,  the  Bhattacharyya  separation  decreases  as 
the  intensity  of  the  source  becomes  weaker.  When  the  source 
intensity  is  much  larger  than  the  background  variance,  the 
Bhattacharyya  distance  is  not  affected.  As  the  range 
increases,  a  point  is  reached  where  the  Bhattacharyya  separ¬ 
ation  falls  off  abruptly.  This  is  the  range  at  which  the 
source  signal  would  be  hidden  by  the  background  and  the 
range  estimate  obtained  would  be  unreliable. 

If  the  range  increment  between  distributions  is 
increased,  the  Bhattacharyya  distance  will  also  increase. 
The  range  increment  that  maintains  a  Bhattacharyya 
separation  of  2.0  is  plotted  in  Figure  16.  For  infinite 
intensity  (no  noise),  range  increments  of  less  than  2  km 
will  maintain  Bhattacharyya  separation  of  2.0  until  18  km. 
For  an  apparent  source  intensity  of  10  W/sr  at  2  km,  range 
distributions  separated  by  2  km  maintain  a  Bhattacharyya 
separation  of  2.0  only  to  4  km.  Beyond  5  km,  a  Bhatta¬ 
charyya  separation  of  2.0  cannot  be  maintained  regardless  of 
the  range  increment.  Because  a  Bhattacharyya  separation  of 
2.0  represents  little  overlap  between  distributions,  the 
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Figure  16.  Range  Resolution  of  ARI  Sources  With  Noise 
Added 


plot  in  Figure  16  can  be  considered  to  be  an  estimate  of  the 
range  resolution  that  can  be  obtained  for  the  given 
intensities . 
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V.  Ranging  Algorithm 

Development 

In  the  previous  chapter,  spectral  bands  were  chosen  so 
that  the  separation  between  range  distributions  in  color 
space  was  maximized.  The  resulting  normal  distributions  can 
be  used  to  determine  decision  regions  into  which  unknown 
sources  will  be  classified. 

All  points  were  fitted  to  an  Nth  degree  polynomial  in  e 
using  a  least  squares  method.  The  order  of  the  polynomial 
was  increased  until  the  coefficient  of  determination  ex¬ 
ceeded  0.98.  Only  the  distributions  representing  ranges 
between  2  and  18  km  were  used.  Above  18  km,  the  extinction 
of  the  spectra  becomes  too  large  to  permit  reliable  results. 
Below  2  km,  the  range  distributions  change  rapidly  and  low 
order  polynomials  do  not  accurately  predict  the  functions. 
All  ranges  below  2  km  will  therefore  be  assigned  to  the  same 
decision  region. 

All  data  points  listed  in  Appendix  B  were  plotted  in 
Figure  13.  In  Figure  17,  only  the  means  of  the  range 
distributions  are  plotted.  A  second  order  polynomial  F(e) 
that  can  be  used  to  predict  the  curve  of  the  means  is  also 
plotted  on  the  same  graph.  From  the  data  obtained  with  the 
ARI  files,  the  function  F(e)  was  found  to  be: 

F(e)  =  0.0310  +  (1.9339)e  -  (1.3125)e2 
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Figure  17.  Mean  of  Each  Range  Distribution 

Each  point  on  this  curve  represents  the  location  of  the  mean 
of  a  normalized  distribution  for  a  particular  range. 

By  adding  a  third  dimension  to  the  same  coordinate 
system,  the  range  for  each  point  of  F(e)  could  be  plotted  on 
a  separate  axis.  The  result  is  a  function  in  e  and  $  that 
represents  the  range  of  each  distribution.  However,  <j>  is  a 
dependent  variable  represented  by  F(e),  and  the  function 
would  not  be  defined  for  all  independent  values  of  e  and  4> . 
The  range  function  actually  depends  on  only  one  independent 
variable  and  can  therefore  be  represented  by  a  function 
R(e).  The  range  of  each  distribution  as  a  function  of  e 
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Figure  18.  Range  of  Each  Distribution 

only  is  plotted  in  Figure  18.  The  fourth  order  polynomial 
R(e)  that  will  predict  range  is  also  plotted,  and  is  given 
by: 

R(e)  =  40.9698  -  (395.7266)e  +  (1815.1719)©2 

-  (4032.4529)e3  +  ( 3400. 9910)e4  (19) 

For  R(«)  to  be  valid  though,  the  value  of  e  must  lie  on  the 
curve  of  the  means. 

When  an  unclassified  source  is  plotted  in  reduced  color 
space,  it  generally  does  not  fall  on  the  curve  of  the  means. 
To  calculate  a  range  for  the  point  using  R(e),  the  point 
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Figure  19.  Slope  of  Decision  Boundary  Between 
Distributions 

must  first  be  projected  onto  F(e).  This  projection  should 
follow  the  decision  boundaries  that  separate  the  distri¬ 
butions.  If  it  is  assumed  that  the  covariance  matrix  for 
each  distribution  changes  slowly,  the  decision  boundaries 
for  points  near  the  curve  of  the  means  are  linear.  The 
slopes  of  these  decision  lines  vary  along  F(e). 

A  method  to  calculate  the  slope  of  the  line  separating 
two  distributions  whose  covariance  matrices  are  equal  is 
given  in  Appendix  A.  Using  this  method,  the  slopes  of 
linear  decision  boundaries  between  each  distribution  derived 
in  Chapter  IV  was  calculated.  The  results  are  plotted  in 
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Figure  20.  Projection  of  Unknown  Source  Along  Decision 
Boundaries 

Figure  19.  For  e  values  representing  the  means  of  range 
distributions  between  2  and  18  km,  the  slope  of  the  decision 
boundaries  can  be  represented  by  the  fourth  order  polynomial 
M(e)  also  plotted  in  Figure  19  and  given  by: 

M(e)  =  5.6837  -  (27.7977)e  +  (142.5650)e2 

-  (412.0620)e3  +  (557.0090)e4  (20) 

If  the  values  of  e'  and  4>*  in  Figure  20  represent  an 
unclassified  source  in  reduced  color  space,  the  slope  of  the 
line  between  the  points  (©’,$’)  and  (e,$)  is  given  by: 


m  =  ($'-<)>)/(©'-©) 


(21) 


m  '* 


s 


where 


(©,$)  =  unknown  projection  of  (©',$')  onto  F(e) 

Because  the  point  (©,$)  lies  on  F(e),  the  value  of  $  can  be 
represented  by  F(e).  In  addition,  the  value  of  m  can  be 
represented  by  M(e).  The  result  is  a  fifth  order  polynomial 
in  ©  given  by 


M(e ) 


d>  *  — F  ( e ) 
e7  -e 


(22) 


Solving  this  equation  for  e  and  substituting  in  values 
obtained  from  the  ARI  files, 

0  =  c0  +  cie  +  C2e2  +  C3e3  +  C4e4  +  C,.©5  (23) 

where 


CQ  =  (-5.6837)©'  -  0.0310  + 

C  =  (27.7977)6'  -  1.9339  +  5.6837 

C2  =  (-142.5650)©'+  1.3125  -  27.7977 

C3  =  (412.0620)©'+  142.5650 

C.  =  (-557.0090)©'-  412.0620 
4 

C$  =  557.0090 


The  solution  to  this  equation  is  a  value  for  ©  that  is  used 
in  the  range  function  R(e).  Because  the  equation  is  a  fifth 
order  polynomial,  there  are  five  possible  roots.  However, 
there  is  only  one  real  root  for  values  of  ©  that  represent 
ranges  between  2  and  18  km.  The  polynomial  was  solved  by 
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calling  a  subroutine  in  an  IMSL  library  (Ref  15). 

Using  this  method,  any  value  of  e'  and  can  be  pro¬ 
jected  onto  the  curve  of  the  means  and  a  range  calculated 
using  the  function  R(e).  Because  the  decision  boundaries 
are  not  actually  linear,  projection  accuracy  decreases  as 
the  distance  from  the  curve  of  the  means  increases. 


Testing 

Although  spectral  regions  are  chosen  that  minimize  the 
effect  of  LOWTRAN  model  variations,  some  inaccuracy  is 
expected  when  atmospheric  conditions  are  different  from  the 
U.S.  standard  atmosphere  that  was  used  in  developing  the 
distributions.  Using  LOWTRAN  model  6  (U.S.  standard), 
spectra  simulating  unknown  ranges  were  created  from  the  ARI 
files.  The  predicted  ranges  are  shown  in  Table  III  for 
several  ranges.  The  range  for  all  sources  fell  within  0.8 
km  of  the  LOWTRAN  simulated  range.  This  implies  that  the 
analysis  is  self-consistent  since  the  1962  U.S.  Standard 
atmospheric  model  was  the  only  one  used  to  generate  the 
distributions . 

The  ARI  files  were  then  multiplied  by  LOWTRAN 
transmittance  from  model  2  (mid  latitude  summer)  and  model 
3  (mid  latitude  winter),  and  the  range  estimates  for  these 
files  are  also  listed  in  Table  III.  Using  model  2,  the 
range  estimates  were  consistently  low;  for  model  3,  the 
estimates  were  consistently  high.  These  are  considered 
systematic  errors  and  could  be  accounted  for  relatively 
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Table  III 

Computed  Ranges  For  ARI  Files 


File 


File 


(a)  LOWTRAN  Model  6  (U.S.  Standard) 


5  km 

8  km 

11  km 

14  km 

mi 

4.8 

7.8 

10.8 

14.2 

16.2 

4.6 

7.7 

11.0 

14.4 

16.6 

4.6 

7.6 

10.9 

14.4 

16.6 

5.4 

8.2 

11.3 

14.8 

16.9 

5.0 

7.9 

11.0 

14.5 

16.6 

5.3 

8.1 

11.0 

14.3 

16.4 

5.3 

8.1 

11.0 

14.3 

16.4 

5.4 

8.3 

11.3 

14.7 

16.7 

5.6 

8.4 

11.4 

14.6 

16.7 

(2)  =  less  than  2  km 
LOWTRAN  Model  2  (Mid  Latitude  Summer) 


2  km 

5  km 

(2) 

4.3 

(2) 

4.5 

(2) 

4.5 

2.0 

4.8 

(2) 

4.4 

2.1 

4.8 

2.1 

4.8 

(2) 

4.8 

2.2 

5.1 

8  km 


11  km 

14  km 

9.7 

9.8 

ill 

9.8 

10.2 

iU  cWimS 

9.9 

12.7 

9.9 

12.6 

9.9 

12.6 

10.2 

12.9 

10.3 

12.9 

17  km 


14.9 

14.5 

14.5 

14.4 

14.6 
14.8 


(c)  LOWTRAN  Model  3  (Mid  Latitude  Winter) 


File 


2  km 


5km  8  km  11km  14  km  17km 


12.6 

12.8 

12.7 

13.1 
12.9 

12.8 

12.7 

13.1 

13.1 


.6.1 

6.3 

6.3 

16.6 

16.4 

16.2 
16.1 

16.5 

16.5 


8.0 

8.4 

8.4 

8.7 

8.4 

8.2 

8.1 

18.4 

18.5 


Figure  21.  Location  of  AFG  Files  in  Color  Space 


easily  by  inputting  an  atmospheric  factor  prior  to  using  the 
algorithm.  In  any  case,  the  range  estimates  are  still  good, 
especially  at  close  ranges. 

The  sixteen  spectra  from  the  AFG  library  representing 
various  infrared  sources  at  unknown  ranges  were  then  tested 
using  the  algorithm.  The  reduced  color  space  plot  of  all 
the  sources  is  shown  in  Figure  21  in  relation  to  the  mean 
curve  F(e).  Except  for  four  files  (#1,  #2,  #12,  and  #13), 
all  points  fall  reasonably  close  to  F(e).  Had  they  not 
fallen  near  the  curve,  the  validity  of  the  algorithm  would 
be  questionable.  This  would  imply  that  the  sources  are 


Table  IV 


r 


Range  Comparison  For  AFG  Files 


File 

Actual 

Calc. 

Error 

Number 

Range 

Range 

(km) 

(km) 

(km) 

rj 


substantially  different  from  the  types  that  generated  the 
original  distributions,  and  new  distributions  should  be 
formed  to  include  the  new  sources.  By  doing  this,  the 
coefficients  of  the  ranging  algorithm  would  change.  In  some 
cases,  it  may  be  necessary  to  change  the  band  selection  and 
recalculate  the  ranging  algorithm  entirely.  The  best  way  to 
determine  the  effectiveness  of  the  algorithm  is  to  test  it 
on  as  many  sources  as  possible. 

The  calculated  range  of  each  AFG  file  is  compared  to 
the  actual  range  in  Table  IV.  Ten  of  the  sixteen  files  have 
actual  ranges  that  are  within  2  km  of  the  calculated  ranges. 
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Twelve  are  within  2.6  km.  Only  four  files  resulted  in 
errors  greater  than  3  km.  Of  these  four  files,  one  is  hot 
metal  (#02),  one  is  a  flare  stack  (#04),  one  is  a  blast 
furnace  (#06),  and  one  is  a  coke  oven  (#07).  There  is  no 
consistency  to  the  source  type.  Two  of  the  four  spectra 
were  recorded  at  sites  in  Michigan,  and  two  were  recorded  in 
Indiana.  Other  sources  in  these  areas  provided  good  range 
estimates . 

It  has  been  shown  with  the  ARI  files  that  atmospheric 
conditions  can  cause  systematic  errors.  Because  the  .  AFG 
files  were  recorded  at  different  locations  and  on  different 
days,  it  is  impossible  to  determine  if  the  errors  are  indeed 
systematic.  Most  of  the  files  have  range  estimates  that  are 
longer  than  the  actual  ranges,  but  the  difference  due  to 
atmospheric  variations  is  questionable. 

It  is  suspected  that  errors  from  subtracting  the  back¬ 
ground  spectra  from  the  composite  spectra  contribute  to  the 
range  error.  It  has  been  noted  (Ref  5)  that  errors  of  257„ 
are  possible,  and  variations  of  that  magnitude  certainly 
would  affect  the  results. 

In  the  reduced  color  space  plot  of  the  spectra  (Figure 
21),  three  of  the  four  points  fall  fairly  close  to  the  curve 
of  the  means.  Only  one  point  (#02)  appears  to  be  far  from 
the  curve.  However,  the  point  for  file  #12  also  plotted  far 
from  the  curve  and  resulted  in  an  accurate  prediction.  None 
of  these  reasons  can  justify  why  the  four  files  resulted  in 


inaccurate  predictions. 

The  cause  of  the  error  can  be  determined  through  a 
careful  analysis  of  the  source  spectra.  There  is  a  sharp 
spike  on  all  source  spectra  near  2200  cm-*.  This  is  caused 
by  the  absorption  from  N^O  at  2224  cm-*  (Ref  16:40).  Be¬ 
cause  the  peak  that  remains  is  almost  completely  in  band  2, 
the  ratio  of  band  3  to  band  2  is  small.  The  ranging  algor¬ 
ithm  interprets  this  as  a  source  at  a  longer  range  than  it 
actually  is. 

This  N^O  spike  did  not  occur  when  setting  up  the  •  dis¬ 
tributions  using  ARI  files  and  LOWTRAN.  The  20  cm-*  wave- 
number  resolution  of  LOWTRAN  averaged  the  spike  out  to  a 
smooth  spectral  curve.  The  distributions  were  formed  by 
comparing  bands  in  which  there  was  no  N^O  spike. 

There  are  two  possible  solutions  to  this  problem.  One 
is  to  use  a  higher  resolution  atmospheric  model  and  reform 
the  distributions  so  the  ranging  algorithm  can  be  corrected. 
The  other  possible  solution  is  to  average  the  AFG  spectrum 
over  a  20  cm-*  wavenumber  bandwidth  and  then  test  the  algor¬ 
ithm  on  the  averaged  spectrum.  Either  of  these  methods 
should  improve  the  accuracy  of  the  estimates. 


VI.  Conclusion 


Summary 

The  atmosphere  is  almost  always  considered  a  hinderance 
to  the  detection  of  electromagnetic  radiation.  However, 
through  a  careful  analysis  of  the  effect  of  the  atmospheric 
absorption  on  an  infrared  signature,  it  is  possible  to 
estimate  the  range  of  an  emitter.  Although  the  accuracy  of 
the  method  is  limited,  the  entire  process  is  passive.  Any 
received  signal  at  the  proper  frequency  (2175-2250  cm--*-) 
will  result  in  an  instantaneous  estimate  of  the  source's 
range. 

Spectral  discrimination  methods  were  presented  that 
allowed  the  analysis  to  be  independent  of  source  intensity. 
Distributions  of  spectral  representations  in  color  space 
were  formed  by  comparing  several  frequency  bands  within  a 
received  signal.  Bhattacharyya  distances  were  calculated  to 
compare  the  separation  between  eac,\  range  distribution,  and 
the  bands  were  adjusted  to  optimize  the  separation. 

Using  two  frequency  bands  at  2175-2200  cm-^  and  2200- 
2225  cm-^,  blackbody  range  distributions  were  separated  by  a 
Bhattacharyya  distance  of  2.0  to  a  range  of  10  km.  When 
actual  sources  were  tested  using  these  bands,  the  Bhatta¬ 
charyya  distance  did  not  reach  2.0  for  any  range.  The  bands 
were  adjusted  to  2200-2225  cm-*-  and  2225-2250  cm-^,  and  the 
Bhattacharyya  distance  reached  2.0,  but  only  through  5  km. 


Three  colors  had  to  be  used  to  improve  this  separation. 

The  optimum  spectral  region  for  comparing  the  IR 
signatures  with  three  colors  was  between  2175  and  2250  cm-*. 
Division  of  this  region  into  three  25  cm-*  wide  bands  pro¬ 
vided  Bhattacharyya  distances  of  2.0  between  range  distri¬ 
butions  separated  by  2  km  out  to  ranges  of  18  km.  The  mean 
and  covariance  matrix  of  each  distribution  was  used  to 
develop  a  ranging  algorithm.  By  using  the  algorithm  on  the 
sources  that  generated  the  distributions,  range  estimates 
within  0.8  km  were  obtained.  This  verified  the  validity  of 
the  method.  If  a  source  of  unknown  range  is  of  the  same 
type  as  the  ones  that  generated  the  ranging  algorithm,  an 
accurate  estimate  of  range  should  be  obtained. 

The  atmospheric  model  used  to  simulate  range  for  the 
sources  was  then  varied.  As  expected,  the  range  errors 
increased,  but  the  predictions  remained  within  2.0  km  of  the 
actual  range  for  87%  of  the  files.  The  algorithm  is  rela¬ 
tively  independent  of  atmoshpheric  variations,  but  a  more 
accurate  prediction  can  be  obtained  if  an  atmospheric  factor 
is  applied  to  the  algorithm. 

In  developing  the  analysis,  it  was  found  that  cool  gas 
sources  that  emitted  little  radiation  beyond  the  atmospheric 
absorption  band  could  not  provide  an  accurate  range  esti¬ 
mate.  The  algorithm  depends  on  sources  having  relatively 
constant  intensities  between  2175-2250  cm-*  so  that  all 
intensity  differences  within  the  comparison  bands  can  be 


interpreted  as  range  effects.  Radiation  intensity  of  cool 
gas  sources  begins  to  decrease  rapidly  below  2200  cm~^,  and 
the  algorithm  incorrectly  interprets  this  as  a  range  effect. 

Only  a  limited  number  of  source  spectra  were  available 
to  test  the  ranging  algorithm.  Of  the  sixteen  spectra, 
twelve  provided  range  estimates  within  2.6  km  of  the  mea¬ 
sured  ranges.  Ten  were  within  1.6  km  of  the  specified 
ranges.  The  errors  were  mainly  caused  by  a  radiation  spike 
at  2224  cm  from  N^O.  This  spike  causes  one  of  the  band 
intensities  to  be  excessively  large,  and  the  ranging  algor¬ 
ithm  assigns  a  much  longer  range  to  the  spectrum  than  it 
should.  Because  of  the  resolution  of  LOWTRAN,  the  spike  was 
smoothed  out  when  the  algorithm  was  developed,  so  it  should 
be  possible  to  correct  the  error  by  using  an  atmopheric 
model  with  a  higher  resolution. 

Background  noise  severely  degraded  the  separation  that 
could  be  obtained  from  range  distributions,  and  the  amount 
of  degradation  depends  on  the  source  intensity.  For  intense 
sources,  range  resolution  of  2  km  is  theoretically  main¬ 
tained  out  to  ranges  of  18  km,  but  this  resolution  rapidly 
decreased  as  the  source  intensity  became  weaker. 

Recommendations 

Many  restrictions  had  to  be  placed  on  the  atmospheric 
conditions  that  generated  the  distributions.  By  relaxing 
some  of  these  restrictions,  a  better  idea  of  the  sensitivity 
to  various  climates  and  weather  conditions  could  be 


obtained.  The  path  between  source  and  sensor  was  also 
limited  to  sea  level  altitudes.  Variations  caused  by 
changing  the  path  profiles  could  also  be  investigated. 

Several  sources  did  not  result  in  an  accurate  range 
prediction.  The  N^O  spike  at  2224  cnT^  causes  a  high 
intensity  in  one  of  the  bands  and  the  ranging  algorithm 
interprets  this  as  a  source  at  a  long  range.  Cool  gas 
sources  also  do  not  result  in  accurate  range  estimates.  It 
may  be  possible  to  sort  these  types  of  source  spectra  into 
different  classes  and  use  different  ranging  algorithms  to 
estimate  the  range  of  sources  in  these  classes. 

Unfortunately,  all  of  the  AFG  files  used  to  test  the 
algorithm  were  not  hot  gas  sources,  and  the  ranges  of  these 
sources  was  limited.  More  source  spectra  should  be  obtained 
to  further  test  the  ranging  algorithm.  Sources  representing 
aircraft  and  missiles  at  ranges  up  to  20  km  would  be  ideal. 
The  algorithm  determined  in  this  study  may  work  very  well  on 
hot  gas  sources,  but  this  could  not  be  verified  using  the 
limited  sources  available. 


Appendix  A 

Linear  Decision  Functions 

When  two  Gaussian  distributions  are  located  in  a  two 
dimensional  coordinate  system,  the  decision  boundary 
separating  the  two  distributions  is  quadratic  (Ref  8:53). 
In  Figure  8,  a  quadratic  decision  boundary  was  constructed 
by  plotting  the  intersection  of  contours  of  equal  proba¬ 
bility  of  the  distribution  functions  for  both  distributions. 
If  the  distributions  have  equal  covariance  matrices,  the 
decision  boundary  becomes  linear.  The  slope  of  the  line 
separating  two  distributions  of  equal  covariance  is  derived 
in  this  appendix. 

To  find  points  that  define  the  intersection  of  iso¬ 
probability  contours,  the  probability  functions  of  the 
distributions  are  set  equal.  The  solution  to  the  resulting 
equation  is  the  set  of  points  that  represent  the  same 
probability  of  falling  in  either  distribution.  The  prob¬ 
ability  function  for  a  Gaussian  distribution  is  given  by 
(Ref  10:35): 

pOC)  =  (2rr)"d/2|C|"1/2exp  [  ^(X-M)TC'1(5(-M)  ]  (10) 

where 

d  =  dimension  of  X 

^  —  (XpX2 . ) 

Fi  =  mean  of  the  distribution 

C  =  covariance  matrix  of  the  distribution 
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I C  t  =  determinant  of  C 

_i 

C  =  inverse  of  C 

In  two  dimensions,  the  covariance  matrix  is  given  by: 


rala2 


ral°2  °2 


where 


o  =  standard  deviation 
Inversion  of  this  matrix  yields 


°1  °2 


l-r  -r 
\-al°2 


When  covariance  matrices  are  equal,  equating  the  probability 
functions  is  equivalent  to  equating  only  the  exponents  of 
the  probability  functions.  The  result  is 


(x1-m1)tc  1{x1-m1)  =  (X2-m2)tc_1(x2-m2)  (26) 


Letting  =  (x^y^),  this  equation  can  be  reduced  to 


(x9-x1  )ai  -  (y9-y1)ra1a9 

y  =  - j  (x)  +  Constant 

(x2-x1)ra1a2  -  (y2-yi)ai 


Obviously,  this  is  the  equation  of  a  line  whose  slope  is 
given  by  the  coefficient  of  x.  Knowing  the  location  of  the 
means  of  two  distributions  and  assuming  equal  covariance, 
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Appendix  B 

Computed  Data  for  ARI  Files 


In  Figure  14,  the  locations  of  the  points  representing 
the  color  space  location  of  all  the  ARI  files  at  several 
LOWTRAN  ranges  were  plotted.  The  tabulated  data  in  this 
appendix  is  the  computed  values  that  generated  that  graph. 

All  values  were  computed  using  the  following  three 
bands  to  compare  the  relative  intensities  of  a  spectrum: 


band 

1: 

2200  - 

2225 

cm-* 

band 

2: 

2225  - 

2250 

cm-* 

band 

3: 

2175  - 

2200 

cm-* 

These  bands  are  not  in  increasing  order  because  band  1  and 
band  2  are  the  same  bands  that  were  used  with  two  colors. 
It  should  be  noted  also  that  the  sequence  of  bands  affects 
the  values  that  will  be  computed  for  e  and  <|>.  When  con¬ 
verting  to  spherical  coordinates  from  rectangular  coor¬ 
dinates,  the  values  of  x,  y,  and  z  are  not  interchangeable. 

The  first  column  of  the  tables  is  self-explanatory. 
The  column  labled  LOWT  is  the  LOWTRAN  atmospheric  model  that 
was  used.  "M6M  stands  for  model  6,  and  "R**"  is  the  range 
in  km.  The  next  three  columns  labled  Band  1,  2,  and  3  are 
the  actual  integrated  intensities  of  each  band  in  units  of 
source  intensity  (W/sr).  The  magnitudes  are  different 
because  some  sources  are  brighter  than  others.  The  columns 
labled  Angle  1  and  Angle  2  are  the  values  of  e  and  $  respec¬ 
tively  that  are  used  in  the  distributions.  Even  though  the 
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Table  V 

ARI  Computed  Data 
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Noise  Variance 


By  A-C  coupling  an  IR  detector's  output,  the  average 
signal  from  a  background  becomes  zero.  The  fluctuations 
from  the  background  are  still  present  though,  and  can  be 
represented  by  a  covariance  matrix.  The  magnitude  of  the 
background  fluctuations  has  been  analyzed  by  Sanderson  (Ref 
11),  and  a  brief  summary  is  presented  in  this  appendix  to 
show  how  the  actual  value  to  use  in  the  covariance  matrix 
was  obtained. 


If  F(x)  represents  the  random  distribution  function 
that  defines  the  fluctuations  in  the  background  intensity  in 
a  two  dimensional  scene,  the  auto-correlation  function  A(p) 
is  given  by 


A  ( p ) 


J 


F(x)  F(x+p)  dx 


(28) 


The  Fourier  transform  of  the  auto-correlation  function  is 
called  the  power  spectrum  and  is  denoted  by: 

W (Tc )  =  J*A(p)  exp( -2  7T  ik* p)  dp  (29) 

A  one-dimensional  power  spectrum  is  usually  called  the 
Wiener  spectrum  (Ref  3:XVII,9). 

The  variance  of  a  noise  signal  can  be  obtained  by 


integrating  the  power  spectrum  (Ref  11): 


Through  an  analysis  of  background  spectra  obtained  from 
Environmental  Research  Institute  of  Michigan  (ERIM),  it  is 
found  that  the  power  spectrum  of  most  backgrounds  varies  as 
(Ref  11): 

W(lc)  =  B/fc  (31) 
where  B  is  a  constant  that  depends  on  the  environmental 
conditions. 

The  integral  of  this  function  is  divergent,  which 
implies  that  power  spectra  by  themselves  cannot  give  an 
indication  of  noise.  The  power  spectrum  must  be  multiplied 
by  sensor  characteristics  before  a  value  for  the  variance 
can  be  obtained. 

By  using  a  sensor  whose  response  is  a  window  function 
as  shown  in  Figure  22 ,  any  slowly  varying  component  of  the 
noise  fluctuations  will  be  averaged  to  zero.  The  only 
fluctuations  that  will  be  sensed  by  the  detector  are  those 
that  have  periods  short  compared  to  the  width  of  the  window. 
The  effect  of  this  type  of  sensor  is  therefore  to  reduce  the 
low  frequencies  that  are  part  of  the  power  spectrum. 
Because  the  detector  samples  over  a  finite  size,  only  an 
average  signal  over  the  width  of  the  window  function  is 
measured.  This  causes  smoothing  out  the  rapidly  varying 
frequencies  in  the  noise  signal.  This  effect  results  in  a 
rolloff  in  the  higher  frequencies  of  the  power  spectrum. 


Figure  22.  Typical  Window  Function 


The  power  spectrum  can  now  be  integrated  to  give  a 
variance.  This  is  done  numerically  by  Sanderson  (Ref  11), 
and  the  actual  value  depends  on  the  specific  window  function 
parameters  and  the  value  of  B  obtained  for  the  power  spec¬ 
trum.  One  of  the  detector  types  described  by  Sanderson  was 
chosen  that  resulted  in  a  variance  of  (1.62)B.  A  value  of  B 
equal  to  1.0  was  also  chosen.  The  detector  type  represents 
one  that  results  in  the  least  amount  of  background  noise 
variance,  and  the  constant  B  represents  the  background  that 
presents  the  most  variance  in  background  noise. 

Because  the  data  obtained  by  ERIM  was  given  in  temper- 
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ature  units,  the  variance  also  resulted  in  temperature 
units.  The  conversion  to  radiance  units  can  be  calculated 
by  using  the  differential  form  of  Planck's  law.  The  result 
is : 

<t£  =  (hci/L)  /  (kT2  )  a2  (32) 

where 

h  =  Planck's  constant 
c  =  speed  of  light 
k  =  Boltzmann's  constant 
V  =  wavenumber  of  radiation 
T  =  temperature  of  background 
L  =  background  radiance 
This  result  is  only  valid  for 

(hci/)/(kT)  »  1  (33) 

In  the  IR  region  of  the  spectrum  with  background  temper¬ 
atures  near  300  K,  this  assumption  is  valid. 

The  conversion  from  temperature  units  gives  an  answer 
in  spectral  radiance  units.  To  convert  this  to  intensity 
units,  the  spectral  radiance  variance  is  multiplied  by  the 
field-of-view  of  the  sensor,  the  square  of  the  distance  to 
the  source,  and  the  bandwidth  over  which  the  intensity  is 
measured : 

<r2  =  [FOV  R2  dv  ]2  a2  (34) 
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The  following  input  parameters  were  used  to  obtain  a  result: 


T  =  300  K  (normal  earth  temperature) 
v  =  2200  cm“^  (near  center  of  bands  used) 

FOV  =  10"^  sr  (typical  sensor) 

dv  =  25  cm“^  (width  of  bands  used) 

The  final  result  is  therefore: 

A 

°  =  3.75  x  10~3  R2  (W/sr)2 

This  value  was  used  for  the  variance  in  all  three  bands 


(35) 
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Appendix  D 

Transforms  of  Band  Variation  to  Angular  Variation 

In  Appendix  C,  a  value  for  the  noise  variance  was 
obtained  for  the  three  bands  used  in  comparing  spectral 
intensities.  Because  the  signal  representations  are 
compared  using  spherical  coordinates,  the  variances  given  in 
Appendix  C  must  be  converted  to  variances  in  terms  of  e  and 
<}>.  The  covariance  matrix  for  noise  could  then  be  added  to 
the  covariance  matrix  of  the  range  distributions  to  get  a 
total  covariance  matrix  for  the  distributions. 

Assuming  small  variations,  the  conversion  can  be 
obtained  by  manipulation  of  the  differentials  of  the 
coordinate  transforms.  Using 


p  =  (x2+y2+z2)^ 

(14a) 

e  =  arctan(y/x) 

(14b) 

o  9  9 

<j>  =  arccos [ z/ (x^+y^+z^ )  ] 

(14c) 

differential  of  the  angles  in  terms  of  x. 

y,  and  z  are 

de  =  (de/dx)dx  +  (de/dy) dy 

(36a) 

=  [-y/ (x2+y2 ) ]dx  +  [x/ (x2+y2 ) ]dy 

(36b) 

Similarly, 

d<J>  =  (d<$>/dx)dx  +  (d$/dy)dy  +  (d<|>/dz)dz 
=  xz(x2+y2 )“^(x2+y^+z2 )”^dx 


(37a) 


dc£  =  [(x^+y^)  ^(x^+y  +2^)  •*■]  [  (xz)dx+(yz)dy+(x^+y^)dz] 

(37c) 

These  equations  can  be  used  to  transform  the  variance 
obtained  in  Appendix  C  into  variances  in  terms  of  e  and  <j>. 
The  magnitude  of  the  source  intensity  (x,  y,  and  z)  must  be 
known  though. 

The  coordinate  transform  equations  allow  calculation  of 
all  the  quantities  given  any  three.  At  any  range,  the  mean 
values  of  e  and  <j>  are  given  in  the  table  in  Appendix  B. 
Therefore  to  completely  specify  all  values  of  x,  y,  and  z, 
only  one  of  the  variables  must  be  given.  For  the  calcula¬ 
tions  in  Chapter  IV,  the  intensity  in  band  1  was  specified. 


y  / 
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